JOURNAL OF RESEARCH of the National Bureau of Standards— D. Radio Propagation 
Vol. 66D, No. 5, September-October 1962 

Scattering From a Conducting Sphere Embedded in a 
Semi-Infinite Dissipative Medium 

Janis Galejs 

Contribution from Applied Research Laboratory, Sylvania Electronic Systems, Waltham 54, Mass. 

(Received March 12, 1962; revised April 18, 1962) 

The conducting sphere is embedded in a semi-infinite dissipative medium and is excited 
by a surface wave or by vertical electric and horizontal magnetic dipoles from the lossless 
half-space. The sphere acts as a combination of horizontal magnetic and electric dipoles. 
Its dipole moment depends on the wave incident from the interface of the two media and on 
the multiple reflections between the sphere and the interface. The transient reflections are 
calculated after establishing the overall harmonic response of the system. 



1. Introduction 

A two-region problem of a condueting sphere in an 
infinite dissipative medium has been considered ])v 
Wait [1951, 1953, 19()()a| for an excitation by a time- 
varying magnetic field, by a coil, which carries an 
oscillating current, or by dipoles. In this note the 
above problem is extended to a condncting sphere 
embedded in a semi-infinite dissipative medium. 

Related topics, applicable to the solution of the 
overall problem, have been discussed in literature. 
Fields of dipoles above and in lossv media [Soin- 
merfeld, 1949; Stratton, 1941; Norton, 1936, 1937; 
Baiios and Wesley, 1953, 1954; Wait, 1961], the 
transient propagation in a lossy jnediuin [Richards, 
1958; Keilson and Row, 1959;* Galejs, 1960; Wait, 
1960c], and scattering by a coiuhicting sphere 
[Stratton, 1941; Wait, 1951, 1953, 19G0a] have been 
considered in detail. Additional references on the 
subject are listed in the above publications. 

This analysis considers the subsurface fields due 
to a surface wave and due to electric and magnetic 
dipole excitations in section 2. The fields scattered 
by the metallic sphere into the lossless medium are 
computed in several steps. In section 3.1 the fields 
scattered by the sphere are computed as in an infi- 
nite dissipative medium, and the sphere is shown to 
act as a combination of a horizontal electric and a 
horizontal magnetic dipole. In a semi-infinite dis- 
sipative medium there will be multiple reflections of 
the fields between the sphere and the interface of the 
two media. These reflections, which are considered 
in section 3.2, provide a correction to the equivalent 
electric and magnetic dipole moments of the sphere. 
The fields of the equivalent magnetic and electric 
dipoles, which penetrate into the lossless medium, 
are computed in section 3.3. After establishing an 
overall response to harmonic excitations in section 
4, the transient response is calculated by Fourier 
transform techniques in section 5. 

2. The Incident Wave 
2.1. A Plane Surface Wave 
The magnetic field near the surface of region 1 is 
assumed to be of the form 



Hf,=Ho{c.)exp(-ikx). (1) 

The field components exhibit a suppressed e'"^ time 
dependence. The definition of coordinates is shown 
in figure 1. In region 2 a compatible solution is 
of the form 

Hll=HM exp (-^fo-72.^■), (2) 

where 
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provided that the displacement curi'cnts of 

2 are negligible (aje<^<^(T). Equations (1) 
approximate tlie variation of the far fields of a dij^ole 
over a hmited range of x. The dimensions of region 

3 are much smaller than the free space wavelength 
\=2T/k and the x variation of the field Ilyo is 
negligible over the sphere of radius a. The field 
at x — O, z = Zo is 

HU{-^o)=HoM exp (-723o). (5) 

The field variation over the volume of th(* sphere 
can be written as 

Hf,iz)=H^Mzo)e-'^^^-^''' ^77j;Uzo)[ l-72(^- .eo)] (6) 
if the radius of the sphere a<5 = V2/(coaio(7). 
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Figure 1. The three region geometry. 
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2.2. Near Fields of a Vertical Electric Dipole 

Fields of dipoles that are located on or above a 
semi-infinite lossy dielectric liave been considered 
by numerous authors [Sonnnerl'eld, 1949; Stratton, 
1941; Norton, 1936 and 1937]. The analysis has 
been usually restricted to fields above or at the 
surface of the lossy medium. Sources located in the 
lossy medium are considered by Banos and Wesley 
[1953 and 1954], while arbitrary coordinates of the 
source and of the observation point are covered in 
the analysis of Wait [1961]. The results of Wait 
[1961] for the horizontal electric dipole will be 
utilized in the present analysis. 

Tlie near fields of a vertical electric dipole (super- 
script ve) can be obtained from the results which are 
available for horizontal electric dipoles (superscript 
he) and from reciprocity relations. It can be 
shown [Wait, 1960b] that the vertical electric field 
of a horizontal electric dipole E^^ is related to the 
radial electric field of a vertical electric dipole 
E^ by 



EV^=E' 



{Idsy cos (t> 



(7) 



where the angle cp is measured from the axis of the 
horizontal electric dipole and where the coordinates 
of source and observation point of E^^ are inter- 
changed relative to coordinates of £'J^ Currents of 
the VE and HE dipoles flow in z^ and x directions, 
where the xyz coordinate system is parallel to the 
x/i/z' coordinate system. The radial electric field 
in the lossy medium of a vertical electric dipole of 
moment Ids, which is located in the lossless medium, 
is obtained from (25) and (30) of Wait [1961] and 
from (7) as 



Idsy2 



LTTCFp 






(8) 



where the coordinates of tlie source are Xs, lis, 2s 
(2,,<0) and wliere 



p'=^{x-x,Y^{y-ysY=^Jx"+y'\ 



(9) 



If the condition |72p'i >>1 is not strictly satisfied, 
(8) becomes inaccurate and the quasi-static approxi- 
mations of Wait [1961] should be used instead. The 
magnetic field in the lossy medium is computed as 






1 ^E. 



(10) 
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Substituting (8) into (1 0) and assuming Zs to be small. 
Ids 



TTie 



27rp 



j,{l+kp')e-''^'-'-^'\ 



(11) 



The magnetic field at x^O nd z=Zo is designated in 
the notation of the preceding section as 



Hn^o)=--H^:;.{zo) 



(12) 



if y,=0 and Xs<C^. Over the volume of the sphere 
H'ylz) exhibits the same variation as Hyi{^) in (6). 



2.3. Near Fields of a Horizontal Magnetic Dipole 

An expression for the subsurface H field of a 
horizontal magnetic dipole or of a vertical loop of 
area da can be derived with the aid of reciprocity 
relations from tlie H field of a liorizontal electric 
dipole. It can be readily shown tliat 



E^;^7'-- 



-iiofio sin(/)'Z/J''(</)=7r) 



(IdaY 
{Ids)' 



(13) 



w^here the angles 0' and are measured from the 
axes of the magnetic and electric dipoles respectively, 
and where the coordinates of source and observation 
point of E'p'f' are interchanged relative to coordinates 
of H^'^''. The superscript hm refers to the horizontal 
nuignetic dipole and the superscript he refers to the 
liorizontal electric dipole. For negligible heights of 
the source Zg, a substitution of H'^ of Wait [1961] in 
(13) and of (13) into (10) results in 






1+ {ikp') + {:ikp'Y]e-'''^'-^-^'' sin0' (14) 



wliere the angle 0' is measured from the dipole axis. 
Equation (14) is valid under the same conditions as 
(8). The magnetic field W'^? will be in the ?/-direc- 
tion at X — if ^,==0 and x^K^O. It will exhibit the 
2:-variatioii of (6) if z^z^. Equality of (14) with 
the corresponding field component of a vertical 
electric (superscript ve) dipole (11) requires for kp^ 
«1 that 

F"'da^mcfy'lp'=P'ds. (15) 

3. The Scattered Field 

3.1. Scattering by the Sphere in an Infinite 
Dissipative Medium 

The incident field exhibits an exponential ^-varia- 
tion in (6), (11), or (14), and it may be approxi- 
mated by a plane wave in the vicinity of the sphere. 
Scattering of a plane wave by a conducting sphere 
in an infinite homogeneous medium has been worked 
out in section 9.25 of Stratton [1941] and can be applied 
to the present problem without formal difficulties. 
The scattered fields are expressed as a sum of dipole 
and higher order multipole modes, but only the 
dipole modes are significant, if the incident field ex- 
hibits only small variations over the volume of the 
sphere (a<5). The n=\ terms (dipole modes) of 
(9.25.4) of Stratton [1941] may be written as 

Er=—^Hy{zo) 7^V>-^ ^1 IH sm e COS0 

72 (72/') " L 72/ J 

(16) 
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H4'-i] 



COS 6 



H;=^H,{zo) ^ '''j^^,^^/''^ «. [i +:^] sill ^ si" ^ (19) 
-^^[i+;^]}«in^ (20) 

— 6i I 1+— 1 cos eX COS (^ (21) 

where ai and &i are given by (9.25.10) and (11) of 
Stratton [1941]. If |72<^|«1 and c7-c>>o', ai and Ji 
are approximated by 

«,=-(M)!ri+5(,,.)^-3^+3(-^Y+ . . .1 

(22) 
?>:=|(72a)'[l-j|(72«)^-^(T3«)+ . . .] (23) 



where y3 = -<Jiwfiac and where the last term of (23) 
is usually negligible. In the limit of a higlily con- 
ducting small sphere these approximations result in 



ai = -6i/2=-(r2a)V3. 



(24) 



The field components proportional to bi are identical 
to those of a ^-directed electric dipole [Stratton, 1941 ; 
Wait, 1960a] of a dipole moment 



Id^ = Q7rE^{Zo)bi(a2+icoe2)/r2 



(25) 



The field components proportional ai are identical to 
those of an 77-directed magnetic dipole [Wait, 1960a] 
of a dipole moment 



Kdr] = QTioiiJLoHy(2o)ailyl. 



(26) 



The electric dipole moment Id^ of (25) can be attrib- 
uted to the incident electric field, while the magnetic 
dipole moment Kdyj of (26) is thought to be caused 
Iby the incident magnetic field [Wait, 1950]. This 
separation of the effects of the electric and magnetic 
fields is particularly useful if the sphere is excited by 
I a non-plane wave. 

The expressions (16) to (21) will apply as long as the 
local fields are not perturbed by surface reflections. 



However, the surface reflections nmst be considered 
when computing the scattered fields in the vicinity of 
the interface between media 1 and 2. 

3.2. Multiple Reflections Between the Sphere and 
the Surface 

The fields in the vicinity of the sphere may be 
perturbed by surface reflections. If the radius of 
the sphere a is small rela tive to th e skin depth of the 
dissipative medium 5=v2/(<^Aioo'), the effect of these 
second order reflections can be shown to be small. 

The £' and //fields of (16) to (21) will be reflected 
downwards from the interface between the media 1 
and 2, but only the field components which leave the 
sphere at an angle d^w arrive back into the vicinity 
of the sphere. At this angle El2=H'^2=0, the E 
field is in the x direction and the H field is in the y 
direction. If |72|»A', the field incident from the 
sphere E^t and Hyi and the field reflected from the 
surface Exr and Hyj. are related at the surface by 



-L/xr — -t^xi 



(27) 
(28) 



The field components Exr{z^ and Hyr(zo), which are 
reflected from the interface into the vicinity of the 
sphere, may be related to the field component 
Exi^o) and Hy{z^ of the plane wave which excited 
the sphere originally. Equations (17), (18), (20), 
(21), (27), and ^(28) result in 

TT(^ x =9exp (-272^0) 



{^'-^^(2-^ 



"(2723o)d (2723o)'J ^^^^ 

The ratios (29) and (30) apply also to amplitudes of 
successively higher order reflections in the vicinity 
of the sphere. When accounting for the higher order 
reflections, the net fields incident on the sphere are 
given by 

■£'x(^o)|Total = £'a:(^o)X; I ^""'/^"x (31) 

//.(^O) |TotaI=-^.(^0) g [§J^] '• (32) 

The eft'ective dipole moment of the sphere is related 
to the total incident field in the same way it was 
related to the field component of the plane wave in 
(25) and (26). These effective dipole moments of 
the sphere are obtained by substituting (31) and (32) 
for Ex{z^) and Z/^(,eo) in (25) and (26). The eftec- 
tive dipole moment can be expressed as 



71 = 



(33) 



640618—62- 
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where 



n=0 



Fe=-QTbJyl 
i^^=67^^co/xoal/72 



(34) 



(35) 



(36) 



Gi={y2ay iy- exp (—272^0) 



{l['+2i;]+(5fe'} <"' 






(38) 
(39) 



The series (33) and (34) are rapidly convergent for a 
<5==V2/|7L even if the exponential factor of (37) is 
of the order of unity. It follows from section 2.1 
that the variation of the incident field over the 
volume of the sphere may be neglected if a is small 
relative to 8. The neglect of multiple reflections 
between the sphere and the surface is therefore 
permissible {\Gi\<C<C^), if the field variation is small 
over the volume of the sphere, as it has been assumed 
in deriving the field components of section 3.1. 

3.3. Fields of the Subsurface Horizontal Electric 
and Magnetic Dipoles 

The fields of a subsurface horizontal electric 
dipole can be obtained from the results of Wait 
[1961]. The time rate of change of the magnetic 
induction is approximated for |72p|>>l by 



BH 



' 27rpV2 



[l+ikp+{ikpy]e-'^p-''2'o cos (t>, (40) 



which is numerically equal to the voltage induced in 
a vertical loop of 1 m^ area. 

The fields B^"^ produced in the lossless medium by 
the subsurface horizontal magnetic dipole can be 
related to the subsurface fields of a vertical electric 
(superscript ve) and a horizontal electric (superscript 
he) dipole that are located in medium 1 . It has been 
shown b}^ Galejs [1961] that 



EI^(2s=Zo,Zr-^Q)=C^Sc^Hmz: = 0,2'r = -o) 



Kd-n 
Idz' 



(41) 



In the above equation the subscripts s and r refer to 
source and receiver coordinates respectively. The 
horizontal magnetic dipole is oriented in the y di- 
rection. The unprimed coordinates are defined as in 
figure 1 ; <P^ and p' coordinates are defined in a cylin- 
drical system, centered in that location of the air- 



conductor interface, where the fields of the subsur- 
face horizontal magnetic dipole should be detected, 
and with the 0' = O direction parallel to the a:;-axis. 
With Hl',{z', = ^, Zr = ^o) given by (11), it follows 
from (41) that 



^^ " 2p7r V^ikpJ 



e-^'^^-^2^ocos(/). (42) 



The electric field E"^"^ is related to H^' by (13) where 
HI' is obtained from (4), (21), (22), (27), and (28) of 
Wait [1961]. This gives 

+ (ikpr]\e-'''-'^'^cosct>' (43) 

The corresponding time rate of change of the mag- 
netic induction is 



Bl 



dp dz 

Kdv 



27rp2 



[1 + {ikp) + (ikpy]e-'^f^-''^'^ cos (t>. (44) 



The combined surface fields are obtained by adding 
(40) and (44). After assuming that the terms pro- 
portional to (ikpy are negligible it follows that 

B4>^-^^^ {l+ikp)e-''^-''^'^ cos «. (45) 
Zirp 

The equivalent magnetic dipole moment {Kd-qY is 
defined as 

io^ixold^ 



{Kdv)'=Kdv' 



72 



(46) 



After applying (33) to (36), (KdrjY becomes 
QTio)p,QHy{zo) 



(KdvY- 



72 



[b,±G^-a,±6i]- (47) 

L n=0 n=0 J 



If it is permissible to ignore the multiple reflections 
and to apply the approximation (24), (Kdrj)^ simpli- 
fies to 

{Kdr])' = 18TiooiJLoHy(zo)ailyl 



- — QiriooiJiQa^Hyizo). 



(48) 



This is three times the magnetic dipole moment in 
(26), which was due to the incident magnetic field 
only. 

4. Response to a Sinusoidal Signal 

The overall response will be calculated for excita- 
tions by a plane surface wave, by vertical electric 
and horizontal magnetic dipoles. For a plane surface 
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wave (superscript;?/) excitation (5), (22), (23), (45), 
and (47) give 



M'-- 



HoicoJaHcofiQ 



■='ii+itpi 



[.?.'••]■ 



-ikp-2y.^z, COS 4 

(49) 



where tlie // vector of the plane wave is in the 
^/-direction and the p and coordinates of the obser- 
vation point are arbitrary. For a vertical electric 
dipole excitation (11), (22), (23), (45), and (47) 
result in 



Bl 



_I{o))dso?io)^Q 
2^^7 



Ln = J 



(50) 



where both the dipole and the observation point are 
located at the same point in the y=0 plane. For a 
horizontal magnetic dipole (or a vertical loop of area 
da) excitation (14), (22), (23), (45), and (47) result in 



Ln = J 



27rp^ 



-2ikp-2y2^Q 



In the above expression Fn is given b}^ 



(51) 



MoTsO^ 



(52) 



where Gi with i=e or m is defined by (37) to (39). 
The summation over n may be approximated by 3 
when neglecting multiple reflections between the 
small sized, higlil3^ conducting spliere and the inter- 
face between media 1 and 2. The attenuation a 
(in dh) of B^(z=2o) relative to -^0(0=0), which is 
due to the exponential factors in (49) to (51) is 
given by 



=20 lo 



giol^ 



-2722o| 



=8.68X2 ^. 

For p<<X, A:p<<l and the three expressions of 
B^ are related as 

, I{co)ds , I{ui)da 
' 2tp^ ' 27rp' ' 



\Bl'\:m:\Br\=HM: 



(54) 



The above J?^ figures are related in the same way 
as the surface fields i/o(w), Hi' and H^"' oi (5), (11), 
and (14). The magnitude of the leading term of 
(51) may be written as 






(-^T> 



(55) 



In a numerical example, the magnitude of Bl"^ 
will be calculated using the following parameters: 
a=2.5 m, 5 = 25 ni, z^=Vl.b m, p^lO^ m, o-^lO'^ 



mho/m (dry ground), and I{oi)da^ 
follows that |i^n=9X10-' V/ml 



:10^ Ami It 



The response of the buried splieres may be simply 
related to the response of the spheres in free space. 
The leading term of the near field response of the 
buried sphere in (49), (50), or (51) is of the form 

15^1=^ exp( — 2£'o/5). The response of the same 
sphere in free space is \Bfj)\=Alm, where 7n=12 for 
identical dipole excitations and m = 6 if tlie former 
surface excitation ^^0(0)) is appHed directly to the 
sphere. The increase in the scattering cross section 
of the sphere in the dissipative medium is therefore 
sufficient to offset the refraction losses for sufficiently 
shallow depths 2^. 

5. Transient Response 

The overall transient response will be calculated 
for an impulse-type plane wave surface excitation 
IIoT5{t), for a step of voltage T^ applied to the elec- 
tric dipole of capacity C and an impulse of current 
Io'T8{f) applied to the magnetic dipole. This results 
with neglected multiple reflections in 

B^{io)) = iccD\ l+io:2)+(i{ic^)-^-^' 

+ qKio^)-'+ . . .] 



-iuih— ■\io} 



(56) 



where 



aiHl 



MoA'MO'c'^' 



A !'/ — 2-\VoO"2;o. 



(57) 



(58) 



The constants D, p, and h have the following values 
for the three types of excitations considered: 
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The transient response is calculated with the aid of 
the Fourier transform pairs (208), (802), (806), and 
(801) of Campbell and Foster [1948]. This results 
in 



'■^^15y') ^ 



where 
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Figure 2. Transient response for p = 0. 
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Figure 3. Transient response for pT^O. 

The response (51) has been plotted in figures 2 and 
3 for several values of the parameters {jpjg) and gi. 
The Vt^Q, g'l^O curves differ insignificantly from 
the p^Q, g'i=0 curves, unless gi>>l. If the terms 
proportional to {pig) and to gi are negligible, this 
response exhibits a minimum of 



5^(r)Uin=^-1.69Z>/^^ 



(63) 



for ^=0.92, it is zero for y=l.^^ and it exhibits a 
maximum 



^0(r)Uax=13.2Z?/^2 



(64) 



for ^max = 4.08. 

The following parameters will be considered by 
the way of a numerical example: a=2.b m, ZQ=l2.b m, 
p=100 m, (7=10~^ mho/m, a-c=10^ mho/m, fi/ /j,o=W, 
and a horizontal magnetic dipole of Ioda=10^ A/ni^. 
It follows that ^1=2.3X10-2 and p/g=0.85. This 
response exhibits a maximum value of B^''?(r)|max 
=430D/^l This is equal to 64X10"^ V/m^ if 
r=10-^sec. 

6. Discussion 

In the preceding calculations, the semi-infinite 
lossy medium was assumed to be excited by a plane 



surface wave or by a vertical electric and horizontal 
magnetic dipole that are located in the interface be- 
tween two media. Modification of the model of 
analysis by buried transmitters and receivers offers 
no further advantages. For horizontal distances p 
much larger than the depth under consideration Zq^ 
the transmitted energy propagates above the lossy 
medium and penetrates vertically downward from 
the interface. A horizontal magnetic dipole located 
in the lossy medium would be coupled inefficiently 
to the surface wave. (A horizontal electric dipole is 
inefficient when its height above the ground is much 
less than X. Its efficiency is only slightly decreased 
for shallow burial depths Zq, where ^o<^=skin 
depth.) Because of the continuit}^ of the tangential 
H field across the interface, a vertical pick-up loop 
would be as sensitive below as above the interface. 
A buried receiver would offer advantages only if it 
would be possible to isolate the wave incident from 
the reflecting object on the interface and the wave 
reflected by the interface. However, even this may 
have merits only in the direct vicinity of the re- 
flecting object (00 >p). 
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